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ABSTRACT 

The rich statistics of galaxy rotations as captured by the velocity function (VF) provides invaluable 
constraints on galactic baryon physics and the nature of dark matter (DM). However, the comparison 
of observed galaxy rotations against cosmological models is prone to subtle caveats that can easily 
lead to misinterpretations. Our analysis reveals full statistical consistency between ~ 5000 galaxy 
rotations, observed in line-of-sight projection, and predictions based on the standard cosmological 
model (ACDM) at the mass-resolution of the Millennium simulation (H I line-based circular velocities 
above ~ 50 km s~^). Explicitly, the Hi linewidths in the Hi Parkes All Sky Survey (HIPASS) are 
found consistent with those in S^-SAX, a post-processed semi-analytic model for the Millennium 
simulation. Previously found anomalies in the VF can be plausibly attributed to (1) the mass-limit of 
the Millennium simulation, (2) confused sources in HIPASS, (3) inaccurate inclination measurements 
for optically faint sources, and (4) the non-detectability of gas-poor early-type galaxies. These issues 
can be bypassed by comparing observations and models using linewidth source counts rather than 
VFs. We investigate if and how well such source counts can constrain the temperature of DM. 



1. INTRODUCTION 

Mass and angular momentum are crucial galaxy prop- 
erties, since their global conservation laws constrain the 
history and future of galaxy evolution (jBullock et al.l 
l2001bl lal). Moreover, measurements of mass and an- 
gular momentum uncover hidden dark matte r and 
potentially constrain its nature (Zavala ~et al.l 120091 : 
[Obreschkow et al. 2013). In recent decades, the mass 
statistics h ave been studied i n detail via the mass func- 
tion (MFjLT^WWtiJp^ the luminosity function 
(LF,|Lovedav et al."2012^7 

and the auto-correlation func- 
tion ([Blake et al. 2011). By contrast, angular momen- 
tum remains a side-topic, normally addressed indirectly 
via the Tully-Fisher relation (TFR, McGaugh 2012) or 
used as a means of recovering the ma ss distribution in 
individual galaxies (|de Blok et al.|[2QQ8t) . Spatial statis- 
tics of angular m omentum and the related circular veloc- 
ity function (VF, i Gonzalez et al.'"2000'; "Desai e t al.ll2004l : 
IZwaan et al.l 12010 : Papastergis et al. 2011) remain rela- 
tively unexplored. This is despite the fact that the VF 
offers a tremendous potential with reg ard to comparin g 
LFs obtained in different wave-bands (jGonzalez et al.l ). 
measuring various mecha nisms of fee dback in the evo- 
lution of galaxies (Sawala et al.l 120121). and constra ining 
the temperature of dark matter (Zavala et al. "2009). 

In fact, measuring a galaxy's rotational velocity is 
challenging, since it requires both a measurement of 
the galaxy inclination, typically drawn from a spatially 
resolved optical image, as well as a measurement of 
the line-of-sight rotational velocity, typically obtained 
from the Doppler-broadening of the 21 cm emission line 
of neutral hydrogen (Hi). Today, only two Hi sur- 
veys offer reasonably large samples to construct VFs, 



the Hi Parkes Ah Sky Survey (HIPASS, iBarnes et all 
120011 ) a nd the ongoing Arecibo Legacy Fast ALFA (AL- 
FALFA, [Gio^nenrit_an 11005^^ They are the largest 
surveys by the cosmic volume and by the number of 
galaxies, r espective ly. The VFs derived from HIPASS 
( Zwaan et all 120101] and the 40%-release of ALFALFA 
( Papastergis et a]Tl2011f ) were both compared against 
theoretical models, including predictions by the S^-SAX- 
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Fig. 1. — (Color online) Counts of Hi linewidths in the largest 
equivalent subsamples of the HIPASS survey and the ACDM-based 
S^-SAX-model. This plot shows that measured and simulated 
linewidths agree on the completeness domain of Mhi ^ 10^ 
and Vc ^ 50 km s"-*^. Error bars represent 67%-measurement un- 
certainties associated with actual measurement noise, cosmic vari- 
ance, and completeness uncertainties. Similarly, the grey shading 
depicts 67%-confidence intervals for the model, associated with cos- 
mic variance. This figure is a simplified version of Fig. (Sj). 
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model (|ObreschkQw et all2QQ9aI ), the only current model 
of frequency-resolved H l-emission lines in a cosmological 
simulation. These comparisons uncovered statistically 
significant differences, some of which could be attributed 
to gas-poor massive early- type galaxies (Zwaan et al.), 
but the physical implications remained unclear. Differ- 
ences in the faint-end of the velocity function (Fig. 9 in 
iPapastergis et all ), near the resolution limit of the S^- 
SAX-model, seemed to hint a possible breakdown of the 
current cosmological model. In a new attempt to under- 
stand and exploit these differences, we successively found 
them to be subtle artifacts of the comparison itself, hence 
motivating a more detailed analysis. 

This paper presents a revised comparison between the 
Hi line profiles in HIPASS and S^-SAX. We deliber- 
ately focus on HIPASS, while reserving a similar anal- 
ysis of the ongoing ALFALFA survey for the future, be- 
cause HIPASS already has optical inclinations available, 
exhibits a detailed completeness function, and contains 
less cosmic variance than 40%-ALFALFA in terms of 
the redshift-distrib ution of the galaxies (s ee Fig. 4a by 
iMartin et al.ll2Ql^ versus Fig. 2 bottom bv I Zwaan et al.l 
[2QQ5I). The HIPASS data is compared against the S^- 
S AX- model in various ways. A key result, worth high- 
lighting early, is the full consistency between the 50- 
percentile Hi linewidth W^q in HIPASS and S^-SAX, 
as illustrated by the counts in Fig. [H In this work we 
compare both apparent Hi linewidths and inclination- 
corrected circular velocities using source counts, as well 
as space density functions. The different aspects uncov- 
ered by these functions are discussed in detail, as well 
as their reliability as statistical estimators. Based on the 
results, we finally conjecture that linewidth counts might 
be a useful tool for measuring the temperature of dark 
matter, and discuss how well HIPASS can, in principle, 
constrain this temperature. 

The manuscript is organized as follows. Section [2] 
first explains the observed dataset (HIPASS with opti- 
cal imaging) and its simulated counterpart (S^-SAX). 
Five statistically independent simulations are generated 
specifically to assess the effects of cosmic variance. The 
observed and simulated datasets are then truncated 
to congruent subsamples suitable for their comparison. 
This comparison is presented in detail in Section [3l In 
Section m the consistency between HIPASS and S^-SAX 
is interpreted and discussed with respect to the TFR and 
alternative models of dark matter. Section [5] summarizes 
the results in a list of key messages. 

2. DATA DESCRIPTION 

2.1. Observed data: HIPASS 

HIPASS is a blind search for Hi emission at declina- 
tions dec < +25° in the velocity range —1, 280 km s~^ < 
cz < 12,700 km s~^, where c is the speed of light and 
z is the redshift. This survey resulted in 5317 identified 
galaxies, gathered in two catal ogs: the HIPASS galaxy 
catalogue fHICAT. lMever et al.|[2QQll Zwaan et al.ll2004 
containing 4315 sources with dec < -h2°, an d its north- 
ern extension (NHICAT, iWong et al.l [20061 ) containme; 
1002 sources with -h2° < dec < +25°. The Hi lines 
of these 5317 sources have been parameterized in vari- 
ous ways. In this work, we will use the luminosity dis- 
tance given in Mpc, the velocity-integrated line fiux 



Sint^ given in Jy km s ^, the corresponding Hi mass 
Mhi = 2.36 • lO^SintDld + ^)-\ given in Mq, the 
peak-fiux density Sp^ given in mJy, and the linewidth 
^50 ('W5'o^'"' in HICAT), given in km s~^ and mea- 
sured at 50% of the peak fiux density. HIPASS uses a 
channel width of 13.2 km s~^, but parameterization was 
carried out after two stages of smoothing (Tukey and 
Hanning), resulting in a full- width-half- max resolution 
o f 26.4 km s" -^ for V^so- 

iDovle et"all C2005) presented optical counterparts for 
HICAT, identified in t he 6,/-band plates o f the Super- 
COSMOS Sky Survey (jHamblv et al.ll200Tl) . To each of 
these galaxies they fitted an ellipse to measure the semi- 
major axis a, the semi-minor axis 6, and the position an- 
gle. There are 3618 sources in HOPCAT with identified 
values a and b. From those values the galaxy inclinations 
i can be estimated using the spheroid assumption. 
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where q = b/a and qo denotes the intrinsic axis ratio, 
here taken to be qp = 0.2 to remain consistent with 
I Zwaan et al.l ([2010). As in the latter work, we here define 
the circular velocity Vc of a galaxy as 
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although the actual asymptotic rotational velocity may 
slightly differ from V^. 

2.2. Simulated data: S^-SAX 

This section summarizes S^-SAX (lObreschkow et al] 
2009a), the first cosmological model of resolved Hl- 
emission lines of galaxies. 

S^-SAX builds on mo del-galaxies generated by a semi- 
analytic model (SAM, iDe Lucia fc Blaizot' '200/ 1). The 
latter relies on the Millennium simulation (Spring el et al.[ 
2005) that tackles the evolution of cold dark matter 
(CDM) in a comoving box measuring (500 /i~^Mpc)^, 
where h is defined via the local Hubble constant Hq = 
100 /ikms~^ Mpc~^. This simulation uses the standard 
cosmological model (ACDM) with parameters h = 0.73, 
Qm 0.25, 0.045, Qa 0.75, and as 0.9. From 

this simulation CDM halos and their merging histories 
are extracted. The SAM then assigns galaxies to the 
centers of these halos using a cooling model and evolves 
global galaxy properties, such as stellar mass, gas mass, 
and morphology according to physical rules allowing for 
feedback from black holes and super novae. The free pa- 
rameters in this SAM were tuned to the locally observed 
color-magnitude distribution, but there is no explicit fit 
to galaxy sizes, rotations, and gas properties. The num- 
ber of model-galaxies at a time of 13.7 • 10^ yr, i. e., 
today, is to about 3 • 10^. Although the cosmological pa- 
rameters of the Millennium simulation are slightly incon- 
sistent with the newest estimates (Komatsu et all[201l[ ). 
the present-day galaxy properties remain nearly unaf- 
fected according to calculations by ?. 

Given the evo l ving model-galaxies of the SAM, 
[Qbreschkow et al] |2009a) assigned refined cold gas prop- 
erties to each galaxy. Their method, sketched out in 
Fig. [21 can be summarized as follows. The scale radius 
of galactic disks is calculated directly from the spin of the 
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Fig. 2. — (Color online) Illustration of our method to model 
H I emission lines. For each model-galaxy, we compute an axially 
symmetric H I density profile Eh i (t) and a circular velocity profile 
y^total^^j based on mass distribution in the halo, disk, and bulge. 
Convolving Eh i (r) with V^°^^^ (r) yields an ideal edge-on emission 
line (panel c, dashed line), which is smoothed with a Gaussian Ker- 
nel of cr = 8 km s"-*^ to account for the velocity dispersion typical 
for the local universe (panel c, solid line). 

dark matter halo. To do so, a variable ratio ^ G [0.5 — 1] 
between the specific angular momentum of baryons and 
dark matter was adopted. This ratio is a function of the 
Hubble- type and the stellar mass, adjusted such that the 
resulting disk scale radii optimally reproduce those of the 
real galaxies inThe HI Nearby Galaxy Survey (THINGS, 
IWalter et al.ll2QQ"8t) . Given the disk scale radius and the 
total cold gas mass from the SAM, the radial H I sur- 
face density Shi(^) is calculated using a pressure-based 
model for the ratio between molecular (H2) and atomic 
(Hi) hydrogen ( Obreschk ow fc Ra wlings 2009), derived 
from the THINGS sample (jLerTv et al. 2008 ). In par- 
allel, circular velocity profiles V^^'^^^r) = (^^^''^"^(r) + 
Ydisk'^I^^^^^Yhu\ge'^^^^^y/2 calculated from the circular 
velocities implied by the gravitational potentials of the 
dark matter halo, the galactic disk, and the central bulge, 
respectively. Convolving Y^'=>^^^(^r) with Shi(^) then re- 
sults in a model for the frequency-resolved Hi emission 
line (dashed line in Fig.[2]3), which, when convolved with 
a Gaussian Kernel for dispersion, turns into a smooth 




Fig. 3. — (Color online) The five colors represent our five virtual 
HIPASS volumes fitted inside the cubic box of the Millennium sim- 
ulation. The box obeys periodic boundary conditions, such that 
the truncated volumes (blue and purple) are in fact simply con- 
nected. The five HIPASS volumes do not overlap and can hence 
be used to estimate the magnitude of cosmic variance. 



profile (solid line). 

Departing from the cubic box of the Millennium sim- 
ulation populated with model-galaxies with resolved 
H I emission lines, Obreschkow et al. (2009b) produced 
a sky- model with apparent extra-galactic Hi emission 
as seen by a fixed observer. To do so, the Cartesian 
coordinates (x, z) of the model-galaxies were mapped 
onto apparent positions, i.e., right-ascension (RA), dec- 
lination (dec ) , and redshift using the method of 
Blaizot et al.' (2005). This method explicitly accounts 
for the fact that galaxies more distant from the observer 
are seen at an earlier stage in their cosmic evolution. 
Along with this mapping, the intrinsic H I luminosities 
are transformed into observable fluxes. Moreover, the 
Hi emission line of each galaxy is corrected for the in- 
clination of the galaxy, respecting, however, the isotropy 
of turbulent /thermal dispersion. The linewidth W50 at 
the 50%-level of the peak flux density is then measured 
from the apparent Hi line of the inclined model-galaxy. 
To allow a clean comparison with observations, we then 
calculate the circular velocity Vc of a model-galaxy via 
eq. ([2j). This is an important step, since Vc can differ from 
the asymptotic value of ^^0*^1^^ _^ up to 30% for 

some galaxies with relatively compact H I distributions. 

For the purpose of this paper we realized five differ- 
ent virtual skies by placing the cosmic volume probed 
by HIPASS five times inside the simulation box of the 
Millennium simulation as shown in Fig. [3l There is no 
overlap between these five sub- volumes, making them (al- 
most) statistically independent. These five virtual sky 
volumes will be used to quantify the effects of cosmic 
variance, that is the random effects of the locally inho- 
mogeneous large scale structure. 

In the aim of comparing the S^-SAX-model against 
HIPASS it is crucial to note that the gridded beam of 
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Fig. 4. — (Color online) Simulated galaxies with centers separated 
by less than the FWHM of the HIPASS beam of 15.5' (left panel) 
and with H I emission lines overlapping in frequency at the 20- 
percentile level (dashed lines, right panel) are regarded as confused. 
For comparison with HIPASS, such confused sources are considered 
as a single sources with an emission line (solid line) obtained by 
co-adding the individual constituents. 

the HIPASS data measures 15.5' at full-width-half-max. 
Using the S^-SAX sky-model, we find that this limited 
spatial resolution implies a non-negligible probability for 
two or more Hi disks to be confused, i.e., to fall inside 
the same beam and simultaneously overlap in frequency. 
This confusion must hence be accounted for when com- 
paring observations against simulations. We do so by 
merging simulated galaxies, whose centroids are sepa- 
rated less than 15.5' and whose Hi lines overlap in fre- 
quency. The common H I mass is then taken as the sum 
of the individual components and VF50 is measured from 
the combined line as shown in Fig. HI We further de- 
fine the circular velocity Vc of the merged object as the 
H I mass- weighted average of the circular velocities of the 
components. This procedure reduces the number of sim- 
ulated sources by about 2%. 

2.3. Sample selection 

Let us now construct subsamples of sources in HIPASS 
and S^-SAX using identical selection criteria. Two types 
of samples will be considered depending on whether opti- 
cal counterparts are required for HIPASS sources. These 
counterparts are needed when considering circular veloc- 
ities Vc^ since those require estimates of the galaxy incli- 
nations. However, masses Mhi and linewidths W^o do 
not require optical data. We shall call the larger galaxy 
sample, in which optical counterparts are irrelevant, the 
"reference sample". A subsample of this reference sam- 
ple, in which all galaxies have optical inclinations and 
thus estimates of T4, is then called the "l^-sample". The 
precise selection criteria of these two samples are listed 
in Tab. [Hand explained in the following. 

2.3.1. Reference sample 

Volume truncation: Since the volume of S^-SAX ex- 
ceeds that of HIPASS, the former must be truncated to 
the field-of-view (FOV) and redshift range of HIPASS, 
i.e., dec < +25° and cz < 12, 700 kms~^. These criteria 
truncate S^-SAX to sub-volumes matching the colored 
regions in Fig. [3l 

Mass limit: The limiting Hi mass, above which S^- 
SAX can be considered complete, is about 10^ /ifo^M© 
(lObreschkow e t al. 2009a), where h^o = i.e., /iro = 

1 if Hq = 70 /ikms~^ Mpc~^, as consistent with current 
observations (iJarosik et al.ll2011[ ). The simulated Hi MF 



Selection criterion 


trims HIPASS 


trims S^-SAX 


Reference sample: 






dec < +25° 


no 


yes 


cz <= 12, 700 kms ^ 


no 


yes 


Mhi > 10^ Vo'M© 


yes 


yes 


^50 > 25 kms-^ 


no 


yes 


C(5p,^int) >0.5 


yes 


yes 


C(6'p, 5'int) ^ R 


no 


yes 


Additional criteria for 


Vc-sample: 




dec < +2° 


yes 


yes 


0.86 > R 


no 


yes 


i > 45° 


yes 


yes 



TABLE 1 

The upper list shows the criteria applied to all data presented in 
this paper. The lower list shows the additional selection criteria 
applied when comparing circular velocities Vc, which require 
optical inclination measurements. Here, R G [0, 1] denotes a 
random number drawn separately for every galaxy and equation. 



drops rapidly below this limiting mass due to the limited 
mass resolution of the Millennium simulation. We must 
therefore limit the comparison between simulated and 
observed sources to the mass range Mhi > 10^ H^^Mq. 
This criterion removes 71 nearby galaxies from HIPASS, 
that is about 1.3% of the 5317 identified sources. 

Limiting linewidth: HIPASS does not resolve sources 
with linew idths sma ller than 25 km s~^ (see Section 3.1 
in Mever et al.ll2004 ) . For correctness we therefore apply 
the selection criterion W^o > 25 km s~^ to S^-SAX, al- 
though this only reduces the number of simulated sources 
by about 0.1%. 

Completeness limit: In HIPASS, real sources are de- 
tected with a probability approximated by 

C(5p,5int) =E[pi(Sp-p2)]E[p3(Sint-p4)], (3) 

where E(x) = max{0, erf(x)}. The parameters are p = 
(0.036,19,0.36,1.1) for dec < +2° (jZwaan e^ll [2QQ1 
and p = (0.02,5,0.14,1) for dec > +2° (|Wong et al.l 
[2006h . The same completeness function must be ap- 
plied to S^-SAX. This is done by drawing a random 
number R uniformly from the interval [0, 1] for ev- 
ery simulated galaxy, and retaining the galaxy only if 
C{Sp, Sint) ^ R- In addition, we must account for the 
fact that the completeness function itself is very uncer- 
tain for CjS y., S\nt) < 0.5 (e.^., Fig. 6 in Zwa an et al.l 



120041 ). As in iZwaan et all (120101 ). we therefore only re- 
tain galaxies with C(5'p, Sint) > 0.5. These completeness 
cuts reduce the total number of simulated sources in the 
HIPASS volume to about 8%, while reducing the number 
of observed sources by 514, i.e., by an additional 9.8% af- 
ter the 1.3% mass cut. This concludes the construction 
of the reference samples. 

2.3.2. Vc-sample 

To compare observed and simulated values of the 
reference sample must be further reduced to a HOPCAT 
equivalent sample, i.e., a subsample with optically mea- 
sured inclinations. 
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Volume truncation: We must exclude the galaxies with 
dec > +2°, for which optical inclinations are not readily 
available in HIPASS. In doing so, the number of observed 
and modeled objects shrinks by roughly 20%. 

HOP CAT completeness: Out of all galaxies in the ref- 
erence sample with dec < +2° only 86% yield optical 
inclinations. Most of the remaining objects lie too close 
to the galactic plane, where the stellar foreground de- 
teriorates extragalactic optical imaging. To account for 
this incompleteness, we reduce the number of simulated 
galaxies to 86% by only retaining the objects satisfying 
0.86 < where R G [0, 1] is a random number. 

Inclination selection: Galaxies with inclinations close 
to face-on exhibit poor inclination measurements, which, 
given their small values of i, result in highly uncertain 
values of when using eq. (j2j). As in Zwaan et al. (2010) 
we therefore only retain galaxies with i > 45°, hence 
reducing the sample sizes by an additional 29%. 

3. COMPARISON BETWEEN HIPASS AND S^-SAX 

Given the identically selected samples of observed and 
simulated galaxies, we shall now compare the statistics 
of the Hi line profiles. This comparison will be carried 
out both at the level of direct source counts (Section [32]) 
and at the level of space density functions (Section 13. 3p . 

3.1. Sample size 

Let us first consider the raw size of the observed and 
simulated samples given in Tab. [2l 

The mean number of sources in the five simulated ref- 
erence samples is about 4926 with a standard deviation 
of 482. This standard deviation is significantly higher 
than the Poisson shot noise of ^ v4926 ~ 70, demon- 
strating the non-negligible effect of large scale structure 
in HIPASS. The number of observed sources in the ref- 
erence sample is clearly consistent with the simulation. 
We therefore expect the normalization of corresponding 
source count statistics and space density functions to be 
consistent between observation and simulation. 

By contrast, the mean number of sources in the simu- 
lated K-samples (about 2013) undershoots the number 
of observed sources by about 339 or 14%. This difference 
is slightly larger than the characteristic value of cosmic 
variance of 267, estimated from the standard deviation 
of the object-numbers in the five simulated l^-samples. 
As argued in Section [321 this moderately significant dif- 
ference between the sizes of the observed and simulated 
V^-samples is at least partially explainable by a small 
fraction of inaccurate inclination measurements in HOP- 
CAT. Those tend to assign high inclinations {i > 45°) 
to objects, which in actual fact have low inclinations 
{i < 45°) and should hence be removed from the ob- 
served K-sample. 

3.2. Source counts 

A refined statistical analysis consists of counting the 
number of galaxies, binned by specific galaxy proper- 
ties. The properties of particular interest are the Hi 
linewidth W^o and the circular velocity Vc defined by 
eq. (j2j). For completeness we also analyze the statistics 
of the H I mass Mh i . The source counts of Mh i and 
1^50 are derived from the reference samples. In turn, the 
source counts of Vc^ which require inclination measure- 
ments, must be performed using the smaller V^-samples. 





reference sample 


l^-sample 


Observation 


4732 


2352 


Simulation 1 


4699 


1839 


Simulation 2 


5268 


2194 


Simulation 3 


5034 


2103 


Simulation 4 


4212 


1642 


Simulation 5 


5416 


2285 



TABLE 2 

Number of sources in each sample. The selection criteria are 
listed in Tab. ^ and the simulated volumes are shown in Fig. [S] 



Figs. [5ti-c show the observed (bars) and simulated 
(lines) counts of Mhi, and Vc^ respectively. The 

grey solid lines correspond to the five individual simula- 
tions, while the black lines represent the geometric means 
of these functions. Variations between the five models 
are due to cosmic variance. The observed source counts 
exhibit several error bars, representing the uncertainties 
described in Tab. [3l Some of these uncertainties are sta- 
tistical, while others are systematic and thus correlated 
across different bins. 

The observed and simulated Mhi counts in Fig. [5^ 
are moderately consistent. Four of the five models and 
the mean model show a slight bump around Mhi ^ 4 • 
IO^Mq. This seems to be a feature of the particular SAM 
chosen here, since it is also present in the 6,7-band LF o f 
the same SAM (see Fig. 8 right of iCroton et"alll2006[ ). 
but absent in other SAMs building on the Millennium 
simulation (e.g. Baugh et al. 2005). 

Fig. [5)3 is the central plot of this paper and extends 
on Fig. [H It demonstrates that the simulated linewidths 
1^50 are fully consistent with the observed ones. We em- 
phasize that this consistency requires that the simulated 
and observed samples are constructed according to iden- 
tical selection criteria (see Tab. [1]). Experimenting with 
different completeness functions C further revealed the 
importance of using the smooth completeness function 
C{Sp,Sint) provided for HICAT and NHICAT. A hard 
sensitivity limit, i.e., C{Sp^ Sint) as a step-function, is 
not sufficient in that it induces variations larger than 
the error bars. Moreover, accounting for the confusion 
of sources turns out to be vital. If instead all individ- 
ual galaxies in the simulated sky were considered distin- 
guishable, then the mean source counts are given by the 
dashed line in Fig. [5b. The difference is most pronounced 
at the largest linewidths of W50 ^ 500 km s~^ (artifact 
'2'). Thus, the largest values of W50 in the observed data 
are mostly due to confused sources, i.e., galaxies within 
the same telescope beam and with H I line profiles over- 
lapping in frequen cy space. In constructing the original 
HICAT dataset (Me ver et a l. 2004) an effort was made to 
flag and separate sources exhibiting confused H I line pro- 
files. About 9% of the sources with W^o > 500 km s~^ 
in the reference-sample have been flagged as confused 
(as opposed to 7% in the whole reference sample). By 
contrast, our modelling revealed that most sources with 
^50 ^ 500km s~^ are confused. This means that it may 
be impossible to identify most instances of confusion by 
relying exclusively on the information in the HIPASS 
data. An example of a confused source is shown in Fig.O 
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Fig. 5. — (Color online) Statistical comparison between HIPASS and S^-SAX. Error bars represent various measurement uncertainties 
described in Tab. [3] Thin grey solid lines represent the five statistically independent simulations, while the thick solid lines represent their 
geometric averages. Dashed lines delineate the same averages if source-confusion is not accounted for. Dotted lines represent only late- type 
galaxies, excluding SO and E-types. Green numbers denote artifacts discussed in Section O they match the numbers in the abstract, (a) 
counts of masses Mhi in the reference samples; (b) counts of linewidths VK50 in the reference samples; (c) counts of circular velocities 
Vc = VK5o/(2sini) in the Vc-samples; (d) Hi MP as derived from HICAT by Zwaan et al. (2005) and predicted using all model-galaxies 
in the Millennium box; (e) space density function of W30 for galaxies of all Hubble-types as derived from HICAT by Zwaan et al. (20Tol) 
and predicted using all model-galaxies; (f) Hi VP for galaxies of all Hubble-types with inclinations i > 45° as derived from HICAT and 
HOPCAT by Zwaan et al. (2010) and predicted using all model-galaxies. Since the observed data points use all Hubble-types, the data in 
panels (e) and (f) must be compared against the solid function. The dotted function (only late-types) nonetheless provides a better fit, 
because the predicted class of gas-poor early- types was simply not detectable by HIPASS (details in section [3?3]) . 
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Fig. 6. — Example of a typical confused source in HICAT. This 
object (HIPASSJ1347-30) was assigned a single width of W^o = 
653.6 km s"-*^ (without confusion flag), the highest value of any 
source with S'int > 50 Jy km s-^. The SuperCOSMOS optical bj- 
band image suggests that the H I emission line is a combination of 
two merging systems. 

The counts of circular velocities Vc are shown in Fig. [5b. 
The models are consistent with the observations for Vc > 
50 km s~^, but drastically differ for smaller velocities (ar- 
tifact '3'). The only major difference between Fig.[5]3 and 
Fig. [5t is the inclination-correction [see eq. (|2])]; there- 
fore the excess of observed sources with < 50 km s~^ 
suggests an issue with their inclinations. A systematic vi- 
sual inspection of the 6j-band images of the SuperCOS- 
MOS Sky Survey used in HOP CAT uncovered that a vast 
majority (> 90%) of the galaxies with Vc < 50 km s~^ 
(about 11% of the 2352 objects in the Vc-sample or 6% of 
ah 4315 galaxies in HICAT/HOPCAT) are problematic. 
They are either too faint or too irregular for an optical 
estimation of the inclination, or they simply exhibit er- 
roneous shape parameterisations. Fig. [7] displays three 
representative examples of the latter case. The ellipses 
in Fig. [71 represent the original parameterization in terms 
of minor axis, major axis, and position angle. The axis 
ratios of these ellipses imply inclinations i > 45° [via 
eq. ([1])]. To the naked eye, however, these three galaxies 
are nearly face-on spiral disks {i < 45°), especially in the 
multi-color image of the source HIPASSJ1200-00, which 
is about two magnitudes deeper than SuperCOSMOS. 
Using the 'correct' inclination for this source rather than 
that suggested by HOP CAT, increases Vc roughly by a 
factor two. Since the correct inclination is then below 
45°, this source would be rejected from the 14-sample 
and thus disappear from Fig. [5b. In conclusion, there 
is a small fraction of incorrect shape identifications in 
HOP CAT, which happens to dominate the low-end of 
the Vc counts. Incidentally, this also explains the asym- 
metric scatter skewed towards low rotational velocities 
in t he HO PCAT-based TFR (upper panels in Fig. 3 in 
[Mever et alllM)! . 

3.3. Space densities 

The source counts presented in the previous section 
depend on the selection criteria of the survey listed 
in Tab. [TJ Survey-independent and thus more funda- 
mental statistical measures are the space density func- 
tions (/) = dN/dV. These functions represent the ab- 
solute number of sources, detected or not, per unit of 
cosmic volume and per unit of galaxy properties, such 
as Mhi (Hi MF) or Vc (VF). Evaluating these func- 
tions from empirical data requires inverting the com- 
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Fig. 7. — (Color online) Three examples of the few galaxies in 
HOPCAT (about 6% of all HICAT/HOPCAT objects) with un- 
certain/inaccurate shape parameterizations. The three greyscale 
images are the 6j-band maps from the SuperCOSMOS Sky Sur- 
vey (Hamblv et al. 2001) used in HOPCAT, while the false-color 
image (i-band in red, r-band in green, ^f-band in blue) shows a cor- 
responding deep image obtained b y the Galaxy And Mass Assem- 
bly (GAM A) survey (Driver et aLl l201in . Yellow ellipses represent 
the fits quoted in HOPCAT; they all overestimate the inclinations 
of the galaxies. [Note that position angles in HOPCAT are given 
anti-clockwise from west rather than north.] 



pleteness function, as well as removing the effects of cos- 
mic variance. This is achieved by the two-dimensional 
stepwis e maximum likelihoo d (2D S WML) n iethod devel- 
oped by Zwaan et al." ('2003') and applied bv lZwaan et al.l 
(2005) and Zwaan et al. (2010) to recover the observed 
space density functions of Mhi, ^50, and shown 
in Figs. [SJi-f. Note that the data shown here include 
all Hubble- types. Figs. [5]i-f also display the simulated 
counterparts (solid lines), obtained simply by binning all 
galaxies contained in the redshift z = box of the Mil- 
lennium simulation. This box is large enough for cosmic 
variance to be neglected. However, the observed space 
density functions still obey the same cosmic variance as 
the respective source counts. Therefore the cosmic vari- 
ance uncertainty is plotted with the observed data, al- 
though we derive its value from the variations between 
the five simulated source counts. 

Fig. [5]i reveals that the simulated and observed 
H I MFs are only marginally consistent in the sense that 
the simulation falls within the error bars for about 50% 
of the data points rather than 67%. The fact that the 
agreement was slightly better in source count statistics of 
Fig. [5^ might indicate a minor artifact in the reconstruc- 
tion of the ob s erved H I MF. For example, as suggested by 
I Zwaan et al.l (|2004i) . the 'true' completeness function C 
exhibits a slight dependence on the shape of the H I line 
profile (single-peaked, double-peaked, flat-top) in addi- 
tion to the main dependence on S'p and 5'int. This small 
higher-order effect could be captured by extending the 
2DSWML method to 3D using C{Sp, S int, shape). 

Figs. [5t and [5f suggest clear inconsistencies between 
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Error name 



Type 



Explanation 



cosmology uncertainty {h) systematic 



completeness uncertainty systematic 



cosmic variance statistical 

shot noise statistical 

direct measurement noise statistical 

distance errors statistical 



Variations of the Hubble parameter h on the interval [0.65,0.75]. 
This interval contains the fiducial value of h = 0.7041q q}4 (7-year 
WMAP+BAO. I Jarosik et al.l[20TTh . as weh as t he recently found h = 
0.743 ±0.021 (Spitzer and Hubble space telescope. iFreedman eralll2012[ ). 
Note that it is important to plot /i-related uncertainties, since observed 
and simulated data scale differently with h; e.g., simulated masses are in 
units of /i~^M0, while observed masses are in units of h~^MQ. 
Approximate 67%-confidence intervals associated with systematic uncer- 
tainties of the completeness function C [see eq. (|3j)] in the range C > 0.5. 
We calculate these uncertainties as AC = 0.4(1 — C)/C, which a pprox - 
imately matches the error bars in Figs. 2 and 6 of IZwaan et al.l ()2004f ). 
Note that variations in the completeness function would alter the simu- 
lated data rather than the observed data. In Fig. [5] these error bars have 
only been plotted on the observed data for graphical convenience. 
67%-confidence intervals associated with cosmic variance, as determined 
from the standard deviation between the five simulated samples. 
Approximate 67%-confidence intervals associated with Poisson shot noise, 
calculated as the square root of the number of sources in the bin. 
Approximate 67%-confidence intervals associated with telescope noise 
and limited frequency resolution (Section [2?!]) . 

67%-confidence intervals associated with errors in the spectroscopic dis- 
tance measurement, assuming average line-of-sight peculiar velocities of 
300 km s~^. This uncertainty is largest from nearby sources and therefore 
largest for low-mass galaxies in HIPASS. 



TABLE 3 

Explanation of the different error bars shown in Fig. [S] The error 'type' refers to whole sample. For example, distance errors due to 
peculiar velocities are systematic for an individual source, but statistical at the level of a sample of sources with random peculiar motions. 
In combining multiple errors into a single error bar, statistical errors are added in quadrature, while systematic errors are added linearly. 



the models and observations. In the small velocity range, 
these inconsistencies (artifact '1') directly relate to the 
mass resolution limit of the Millennium simulation. This 
limit implies a significant incompleteness of simulated 
objects with W^o ^ 80 km s~^ and Vc ^ 50 km s~^ (and 
Mhi < lO^Mo to the left of Fig.EJi). In turn, this mass- 
limit is probably linked to the spurious bumps around 
^50 ^ 120 km s-^ and K ^ 70 km s-^ 




Fig. 8. — (Color online) Predicted space density p(Mhi, V^) per 
pixel of size A log^^g ^Hl = ^ logio = 0.1, colored according to 
the average galaxy type in each pixel. Curved lines represent iso- 
lines of the number n of expected detections per pixel; the diagonal 
shading (n < 1) highlights the blind zone of HIPASS. 



A more subtle feature in Figs. [5^ and [5f are the sig- 
nificant deviations at W^o > 500 km s~^ and Vc ^ 
200 km s~^ (artifact '4'). Those deviations are absent 
in the corresponding source counts of Figs.[5]3 and [5b. A 
systematic investigation of the simulated galaxies in this 
high-velocity regime reveals them to be dominated by 
early- type galaxies of numerical Hubble- type T < (E, 
SO) hosting low-mass, but fast-rotating H I disks. Ex- 
cluding those objects from the simulation modifies the 
predicted functions in Figs. [5^ and [Sf to the dot-dashed 
lines, which are in much better agreement with the ob- 
served data, as already noted by'Zwaan et alj (|2010[ ). In 
other words, the model predicts that the high-end of the 
VF is dominated by gas-poor early- type galaxies, but 
it also predicts that HIPASS is unlikely to detect these 
galaxies; hence the consistent source counts. To show 
this explicitly, let us calculate the maximal comoving dis- 
tance I^max (in Mpc) out to which a galaxy {Mhi,K} 
can be detected in the sense that the completeness func- 
tion C drops to 50% at that distance. Substituting 
Sp for 10^5'intK~^ (approximation for i > 45°) and 
5'int for 4.2 • 10~^MHil^max (approximation for z <C 1), 
C(5'p,5'int) = 0.5 [using eq. @] numerically solves to 
^max - 7 • 10-^Mhi exp(-0.4F,0-34) fo^ HICAT and 
NHICAT. The cosmic volume Knax (in Mpc^), in which 
HIPASS can detect a galaxy specified by {Mhi, K} then 
becomes Knax ^ 0.63 • (47r/3)I)^g^x ^ where 0.63 is the 
sky-coverage of HIPASS, i.e., 

V;.ax(MHi, K) ~ 5 • 10-^M^{^ exp {-O.QV°-^^) . (4) 
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On the other hand, the S^-SAX model allows us to pre- 
dict the space-density p(Mhi, K) of a source {Mhi, K}, 
defined as the average number of sources per Mpc^ within 
a pixel {log^LO^Hi i ^/2,log;L0^ ^ (here using 

A = 0.1). The product 

n(MHi, K) = Kiax(MHi, K)p(Mhi, (5) 

then approximates the predicted number of HIPASS de- 
tections per pixel in the {Mhi, K}-plane. Fig. [8] dis- 
plays p(Mhi, Vc) colored by galaxy type with isolines of 
^(Mhi,K)- The region n(MHi,K) < 1 contains less 
than one detection per pixel and thus represents a 'blind 
zone' of HIPASS. This blind zone contains the gas-poor 
(Mhi < 10^ Mo), fast-rotating (K > 200 km s-^) early- 
type galaxies predicted by the model. Since HIPASS is 
very insensitive to these galaxies, it is simply unable to 
recover the predicted high-end of the VF. Surveys deeper 
than HIPASS are needed to verify whether the predicted 
amount of massive gas-poor early- type galaxies is correct. 
For now, it seems safe to conclude that the HIPASS VF 
approximates the VF of late- types, even if no Hubble- 
type cut is applied to the dataset. On a side-note, the 
deeper ALFALFA survey does indeed find significant dif- 
ferences in the hig^h- velocity end of th e velocity function 
(e.g. Fig. 4 in IPapastergis et al.l l2011[ ). 

In principle, the artifacts '2' and '3' of Figs.jSb and[5t 
are still present in Figs.jS^ and[5]F, but they are occluded 
by the even stronger artifacts '1' and '4'. This shows that 
the comparison between models and observations is less 
prone to spurious artifacts, when performed using source 
counts. Furthermore, within the source counts, W^q is a 
less problematic quantity than due to artifact '3'. 

4. DISCUSSION 

This section discusses the physical implications of the 
excellent consistency between observed and simulated H I 
linewidths, as well as potential applications. 



4.1. Interpretation of the consistency 0/VF50 

What does the consistency between the observed and 
modeled Wso-counts (Fig.[5]3) tell us? Does it strengthen 




^50 [kms-^] 



Fig. 9. — (Color online) Analogous plot to Fig.[Tl but with addi- 
tional lines for two alternative dark matter models assuming finite 
particles masses of 1 keVc~^ and 0.5 keVc~^. The error bars sum 
up all the statistical and systematic uncertainties considered in this 
work (see Fig. [Sj and the grey shading denotes the standard devi- 
ation of five independent simulated reference samples for CDM. 



the case of the ACDM model or does it merely manifest 
the empirical tuning of the free parameters in the galaxy- 
model? There is, as argued here, a bit of both. The 
local galaxy stellar MF in the model has been adjusted 
indirectly by tuning the feedback from star formation 
and black holes on the interstellar medium to reproduce 
the o bserved 6j-band and i^-band LFs (jCroton et al.l 
l2006f ). Moreover, the radii of galaxies match the lo- 
cally observed mean stella r mass-to-scale radius relation 
(Obreschkow e t al']l2009af ). One might therefore expect 
the galaxy rotations, which depend roughly on mass and 
radius, to align with local observations. In this argument, 
it should nonetheless be emphasized that the free model- 
parameters (feedback coefficients and the spin ratio of 
baryonic matter to dark matter) have only been varied 
within the restricted ranges consistent with current high- 
resolution observations and high-resolution simulations. 
Therefore, we can at least conclude that the consistency 
of the 1^50-counts in Fig. [SJd confirms ACDM within the 
current uncertainties of galaxy-modelling. 

Moreover, it is worth emphasizing that the relation be- 
tween stellar mass and scale radius is subject to very 
large scatter , both observa tionally and in the model 
(Fig. 2 in lObreschkow et aDl2009af ). Therefore, even if 
the mean relation between stellar mass and scale radius 
is fixed to observations, this merely corresponds to an 
overall normalization of the VF and the corresponding 
Wso-counts. The details of these functions depend on the 
shape of the multi-dimensional probability-distribution 
of halo mass, stellar mass and disk scale radius. This 
shape has not been constrained by empirical fits. Instead, 
it depends directly on the masses, spins, and merging 
histories of the dark halos in the Millennium simulation. 
This argument increases the support of ACDM. 

4.2. Constraints on the dark matter type 

Quantifying the degree to which the Wso-counts sup- 
port ACDM is of course a more delicate affair. For ex- 
ample, what is the actual range of allowed dark matter 
particle masses moM, assumed infinite in CDM but finite 
in Warm Dark Matter (WDM) models? Answering this 
question would require a large array of different WDM 
models, similar to the Millennium simulation, equipped 
with SAMs, where all the uncertainties associated with 
every free parameter are tackled down to the I^so-counts. 
The mammoth numerical requirements of this task lie at 
the edge of current super-computing capacities. 

Here, we limit the analysis to a first order approx- 
imation of the variation of the VFso-counts as func- 
tion of moM, keeping the free parameters of the 
galaxy- model fixed to their best values in ACDM. 
This approximation is obtained by rescaling the num- 
ber density of each galaxy in S^-SAX, one-by-one, by 
^WDM(Mhaio)/^CDM(Mhaio), where Mhalo is the mass of 
the halo containing the galaxy and (/>wDM(^haio) and 
(/>CDM (Mhalo) are the local halo MFs of WDM and CDM 
halos, respectively. These MFs are modeled analytically 
by evol ving the initial density field using the formula- 
tion of lSheth fc TormenI ([TOOOh . WDM models are ob- 
tained by subjecting the initial CDM pow e r spe ctrum 
to a transfer function following 'Bode e tld] (|2001f ) . For 
consistency, these calculations were performed using the 
cosmological parameters of the Millennium simulation. 

Fig. [9] shows the l^so-counts for CDM and two WDM 
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Fig. 10. — (Color online) Relationships between edge-on 
H I linewidths V20 and two mass tracers for local galaxies. A repre- 
sentative sample of 10^ simulated galaxies is shown as black points 
(spiral galaxies) and open circles (elliptical galaxies). Solid lines 
are uniformly weighted power-law fits to the simulated spiral galax- 
ies; in Fig.llOb this fit uses only galaxies with Mgtars > IO^/i^q^Mq. 
The blue dots and dashed lines are observational data and power- 
law fits from McGaugh et al.. (,2000. ) and references therein. 

scenarios with particle masses moM = 1 keVc~^ and 
^DM = 0.5 keVc~^, respectively. Although the observed 
VI/50-counts are only marginally consistent with moM = 
1 keVc~^ and inconsistent with moM = 0.5 keVc~^, 
those WDM cosmologies need not to be incompatible 
with the observed I^so-counts. In fact, we cannot ex- 
clude that varying the free parameters of the SAM within 
the currently allowed ranges can bring the WDM models 
in line with the observed data. However, Fig. [9] conveys 
that if all free parameters in the galaxy-model can be re- 
placed by independently determined precise values, then 
the Wso-counts from HIPASS can indeed discriminate 
between CDM and WDM with 1 keVc~^ particles. 

4.3. Tully- Fisher relation 

So far, we have shown that the H I masses and circular 
velocities of the galaxies in the SAM (as modeled via S^- 
SAX) are consistent with observations; and Croton et alj 
(j2006) showed that the stellar masses are consistent with 
local observations as well. However, the fact that cir- 
cular velocities and masses are independently consistent 
with observations does not, in fact, imply that their two- 
dimensional distribution is correct, too. Therefore, we 
shall finally discuss the two-dimensional distribution of 
circular velocities and baryon masses, i.e., the baryonic 
TFR. To remain consistent with observational standards 
the circular velocity is here approximated as ^20, de- 
fined as half the apparent H I linewidth 1^20 (measured 
at the 20% peak flux level), corrected for inc linations. 
The observational data is drawn from McGaugh et alj 
([200 0) and corrected for h = 0.73. These data include 
galaxy types from dwarfs to giant spirals, whose values 
of V20 have been recovered from H I line measurements, 
corrected for inclinations drawn from optical imaging. 
Only inclinations above 45° were retained to restrict 
the uncertainties of sini. The comparison of these data 
against S^-SAX in Fig. [TOl reveals a good consistency, al- 
though the observational scatter is 50% larger than that 



of S^-SAX. This difference is explainable by measure- 
ment uncertainties, especially regarding the inclination 
corrections in th e low-mass end of Fig. [TOk according to 
iMcGaugh et al.l . Additionally, the S^-SAX-model prob- 
ably underestimates the scatter in V20 by ignoring the 
detailed substructure of H I, such as turbulent mixing in 
mergers, high- velocity clouds, warps, and gas-rich satel- 
lites. 

Unlike the baryonic TFR (Fig. [TOb). the stellar mass 
TFR (Fig. [TOk ) clearly departs from a power-law rela- 
tion for gal axies with ^20 < 200 km s~^. As emphasized 
before (e.g. iMcGaugh et al.l ). this reflects the trend for 
high gas-fractions in low-mass galaxies and conflrms that 
the TFR is fundamentally a relation between circular ve- 
locity and total mass, which is a function of the baryon 
mass (Papaster gis et al. 2012). 

5. CONCLUSION 

This paper presented a detailed comparison between 
the H I lines is HIPASS and those in S^-SAX, a cosmo- 
logical model of galaxies with resolved H I lines. The 
results can be condensed into a list of key messages. 

1. The H I linewidths of the S^-S AX- model are consis- 
tent with those measured from HIPASS (Fig. [SJd). 
Hence, observed Hi linewidths are consistent with 
ACDM at the resolution of the Millennium simula- 
tion (Mm > 10^ M0, Vc > 50 km s"^^ within cur- 
rent galaxy formation models. This does not con- 
tradict a possible breakdown o f ACDM at smaller 
masses (e.g. Zaval a et al.ll2009f ). 

2. Galaxies with Vc < 50 km s~^ tend to be opti- 
cally faint or irregular, thus suffering from large 
inclination uncertainties. To use these objects for 
physical applications, it is better compare simula- 
tions against apparent widths W30 rather than the 
inclination-corrected Vc values. 

3. The model predicts that gas-poor early- type galax- 
ies dominate the high-end of the VF. Yet the model 
also predicts that HIPASS is very insensitive to 
these galaxies because of their small Mhi, large 
W^o (hence higher noise), and low space-density. 
To test whether gas-poor early-type galaxies re- 
ally dominate the high-end of the VF deeper sur- 
veys are needed, but is seems safe to conclude that 
the HIPASS VF obtained using all observed galaxy 
types remains a VF of late- type galaxies. 

4. Most sources with W50 > 500 km s"^ in HIPASS 
are found to be confused; hence confusion must be 
corrected in the high-end of the VF. This flnding 
also applies to ALFALFA, because the ~ 4 times 
higher spatial resolution of the Arecibo beam is 
nearly compensated by the mean redshift being ^ 3 
times higher. 

5. In general, W^o counts are the most reliable statis- 
tics of galaxy rotations^ since they can explicitly ac- 
count for source confusion and complex complete- 
ness functions, and since they are not affected by 
inclinations. On the downside, W50 counts are less 
sensitive to cosmological parameters than velocity 
functions, since each value of W^o mixes galaxies of 
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different masses seen at different inclinations. How- 
ever, the W50 counts of HIPASS are nonetheless 
sensitive to the temperature of dark matter. 

6. In fact, if all free parameters in SAMs can be elim- 
inated or at least constrained independently, the 
l^so-counts derived from HIPASS can verify CDM 
against WDM with 1 keVc~^ particles. 

These cosmological tests and prospects promise to 
become particularly fruitful when applied to future 
Hi surveys, such as the full ALFALFA survey and 
ultimately the ASKAP HI All-Sky Survey (WALLABY) 
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